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Abstract. Spallative nucleosynthesis is thought to be the 
only process capable of producing significant amount of 
Beryllium (Be) in the universe. Therefore, both energetic 
particles (EPs) and nuclei to be spalled (most efficiently 
C, N and O nuclei in this case) are required, which in- 
dicates that supernovae (SNe) may be directly involved 
in the synthesis of the Be nuclei observed in the halo 
stars of the Galaxy. We apply current knowledge relat- 
ing to supernova remnant (SNR) evolution and particle 
shock acceleration to calculate the total Be yield associ- 
ated with a SN explosion in the interstellar medium, fo- 
cusing on the first stages of Galactic chemical evolution 
(i.e. when metallicity < Z < 10~ 2 ). We show that dy- 
namical aspects must be taken into account carefully, and 
present analytical calculations of the spallation reactions 
induced by the EPs accelerated at both the forward and 
the reverse shocks following the SN explosion. Our results 
show that the production of Be in the early Galaxy is still 
poorly understood, and probably implies either selective 
acceleration processes (greatly favouring CNO accelera- 
tion), reconsideration of the observational data (notably 
the O vs Fe correlation), or even new energy sources. 

Key words: Acceleration of particles; Nuclear reactions, 
nucleosynthesis; ISM: supernova remnants; Galaxy: abun- 
dances 



1. Introduction 

There has recently been considerable interest in the Galac- 
tic evolution of the abundances of the light elements Li, Be 
and B (Feltzing & Gustafsson 1994; Reeves 1994; Casse et 
al. 1995; Fields et al. 1994,1995; Ramaty et al. 1996,1997; 
Vangioni-Flam et al. 1998). The Be abundance is partic- 
ularly interesting because this element is thought to be 
produced exclusively by spallation reactions involving col- 
lisions between nuclei of the CNO group of elements and 
protons or alpha particles at energies greater than about 
30MeV per nucleon (MeV/n). Thus the evolution of the 
Be abundance contains information about the particle ac- 
celeration and cosmic ray history of the Galaxy. 



The evolution of the Be abundance, and indeed the 
evolution of all elemental abundances, has to be deduced 
from observations of the fossil abundances preserved in 
the oldest halo stars. Advances in spectroscopy over the 
last decade have greatly improved the quality of the data 
available (Duncan et al. 1992,1997; Edvardsson et al. 1994; 
Gilmore et al. 1992; Kiselman & Carlsson 1996; Molaro et 
al. 1997; Ryan et al. 1994) and the main result is easily 
summarized : in old halo stars of low metallicity, the ratio 
of the Be abundance to the Iron (Fe) abundance appears 
constant, that is to say the Be abundance rises linearly 
with the Fe abundance. 

This has been a surprising result. Naively one had ex- 
pected that, because Be is a secondary product produced 
from the primary CNO nuclei, its abundance should vary 
quadratically as a function of the primary abundances at 
low metallicities. Indeed, considering that the cosmic rays 
(CRs) responsible for the Be production are somehow re- 
lated to the explosion of supernovae (SNe) in the Galaxy, 
it is natural to assume that their flux is proportional to the 
SN rate, dAsN/di. Now since the number of CNO nuclei 
present in the Galaxy at time t is proportional to the to- 
tal number of SN having already exploded, Ns^(t), the Be 
production rate has to be proportional to Ns^dNs^/dt. 
Therefore, the integrated amount of Be grows as Ag N , that 
is quadratically with respect to the ambient metallicity 
(C,N,0 or Fe, assumed to be more or less proportional to 
one another). 

The above reasoning, however, relies on two basic as- 
sumptions that need not be fulfilled : i) the CRs recently 
accelerated interact with all the CNO nuclei already pro- 
duced and dispersed in the entire Galaxy and ii) the CRs 
are made of the ambient material, dominated by H and 
He nuclei. Instead, it might be i) that the proton rich 
CRs recently accelerated interact predominantly with the 
freshly synthesized CNO nuclei near the explosion site and 
ii) that a significant fraction of the CNO rich SN ejecta 
are also accelerated. In both linear growth of the 

Be abundance with respect to Fe or O would arise very 
naturally, since the number of Be-producing spallation re- 
actions induced by each individual supernova would be 
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directly linked to its local, individual CNO supply, in- 
dependently of the accumulated amount of CNO in the 
Galaxy. 

In fact, as emphasised by Ramaty et al. (1997), the 
simplest explanation of the observational data is to as- 
sume that each core-collapse supernova produces on aver- 
age 0.1 M Q of Fe, one to few M© of the CNO elements 
and 2.8 x 1CP 8 M© of Be, with no metallicity dependence. 
Clearly if this is the case and the production of Be is di- 
rectly linked to that of the main primary elements, the 
observed linear relation between Be and Fe will be repro- 
duced whatever the complications of infall, mixing and 
outflow required by the Galactic evolution models. On the 
other hand, although the simplest explanation of the data 
is clearly to suppose a primary behaviour for the Be pro- 
duction, it is possible that this could be an artifact of the 
evolutionary models (as argued, e.g., by Casuso & Beck- 
man 1997). 

Most work in this area has attempted to deduce in- 
formation about cosmic ray (or other accelerated particle 
populations) in the early galaxy by working backwards 
from the abundance observations. While perfectly legit- 
imate, our feeling is that the observational errors and 
the uncertainties relating to Galactic evolution in general 
make this a very difficult task. We have chosen to approach 
the problem from the other direction and ask what cur- 
rently favoured models for particle acceleration in super- 
nova remnants (SNRs) imply for light element production. 
This is in the general spirit of recent calculations of the ir°- 
decay gamma-ray luminosity of SNRs (Drury et al. 1994) 
and the detailed chemical composition of SNR shock accel- 
erated particles (Ellison, Drury & Meyer, 1997) where we 
look for potentially observable consequences of theoretical 
models for cosmic ray production in SNRs. 

Interestingly enough, the study of particle accelera- 
tion in SNRs suggests that both alternatives to the naive 
scenario mentioned above do occur in practice, as demon- 
strated qualitatively in Sect. ||. The first of these alter- 
natives, namely the local interaction of newly accelerated 
cosmic rays in the vicinity of SN explosion sites, has al- 
ready been called upon by Feltzing & Gustaffson (1994), 
as well as the second, the acceleration of enriched ejecta 
through a SN reverse shock, by Ramaty et al. (1997). How- 
ever, no careful calculations have yet been done, taking the 
dynamics of the process into account, notably the dilution 
of SN ejecta and the adiabatic losses. Yet we show below 
that they have a crucial influence on the total amount 
of Be produced, and that a time-dependent treatment is 
required. Indeed, the evolution of a SNR is essentially a 
dynamical problem in which the acceleration rate as well 
as the chemical composition inside the remnant are func- 
tions of time. The results of the full calculation of both 
processes and the discussion of their implications for the 
chemical evolution of the Galaxy will be found in an as- 
sociated paper (Parizot & Drury 1999). Here we present 
simple analytical calculations which provide an accurate 



understanding of the dynamics of light element produc- 
tion in SNRs and elucidates the role and influence of the 
different parameters, notably the ambient density. 

Although Li and B are also produced in the processes 
under study, we shall choose here Be as our 'typical' light 
element, because nuclear spallation of CNO is thought to 
be its only production mechanism, while Li is also (and 
actually mainly) produced through a + a reactions, 7 Li 
may be produced partly in AGB stars (Abia et al. 1993), 
and 11 B neutrino spallation may be important as seems 
to be required by chemical evolution analysis (Vangioni- 
Flam et al. 1996). In order to compare our results with 
the observations, we simply note that, as emphasized in 
Ramaty et al. (1997), the data relating to the Galactic Be 
evolution as a function of [Fe/H] indicate that ~ 1.6 10 -6 
nuclei of Be must be produced in the early Galaxy for each 
Fe nucleus. Therefore, if Be production is indeed induced, 
directly or indirectly, by SNe explosions, and since the av- 
erage SN yield in Fe is thought to be ~ 0.11 M©, each su- 
pernova must lead to an average production of ~ 3.8 10 48 
nuclei (or ~ 2.8 10~ 8 M©) of Be, with an uncertainty of 
about a factor of 2 (Ramaty et al. 1997). We adopt this 
value as the 'standard needed number' of Be per super- 
nova explosion. To state this again in a different way, for 
an average SN yield in CNO of, say, ~ 1 M©, the required 
spallation rate per CNO atom is ~ 3 10~ 8 . 

2. Particle acceleration in SNRs 

It is generally believed that cosmic ray production in SNRs 
occurs through the process of diffusive shock acceleration 
operating at the strong shock waves generated by the in- 
teraction between the ejecta from the supernova explosion 
and the surrounding medium. Significant effort has been 
put into developing dynamical models of SNR evolution 
which incorporate, at varying levels of detail, this basic 
acceleration and injection process (one of the major ad- 
vantages of shock acceleration is that it does not require 
a separate injection process). Qualitatively the main fea- 
tures can be crudely summarised as follows. 

In a core collapse SN the collapse releases roughly 
the gravitational binding energy of a neutron star, some 
10 53 erg, but most of this is radiated away in neutrinos. 
About Esn = 10 51 erg is transferred, by processes which 
are still somewhat obscure, to the outer layers of the pro- 
genitor star which are then ejected at velocities of a few 
percent of the speed of light. Initially the explosion en- 
ergy is almost entirely in the form of kinetic energy of 
these fast-moving ejecta. As the ejecta interact with the 
surrounding circumstellar and interstellar material they 
drive a strong shock ahead into the surrounding medium. 
The region of very hot high pressure shocked material be- 
hind this forward shock also drives a weaker shock back- 
wards into the ejecta giving rise to a characteristic forward 
reverse shock pair separated by a rather unstable contact 
discontinuity. 
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This initial phase of the remnant evolution lasts until 
the amount of ambient matter swept up by the remnant is 
roughly equal to the original ejecta mass. At this so-called 
sweep-up time, isw, the energy flux through the shocks 
is at its highest, the expansion of the remnant begins to 
slow down, and a significant part of the explosion energy 
has been converted from kinetic energy associated with 
the bulk expansion to thermal (and non-thermal) energy 
associated with microscopic degrees of freedom of the sys- 
tem. The remnant now enters the second, and main, phase 
of its evolution in which there is rough equipartition be- 
tween the microscopic and macroscopic energy densities. 
The evolution in this phase is approximately self-similar 
and resembles the exact solution obtained by Sedov for a 
strong point explosion in a cold gas. 

It is important to realise that the approximate equality 
of the energy associated with the macroscopic and micro- 
scopic degrees of freedom in the Sedov-like phase is not 
a static equilibrium but is generated dynamically by two 
competing processes. As long as the remnant is compact 
the energy density, and thus pressure, of the microscopic 
degrees of freedom is very much greater than that of the 
external medium. This strong pressure gradient drives an 
expansion of the remnant which adiabatically reduces the 
microscopic degrees of freedom of the medium inside the 
remnant and converts the energy back into bulk kinetic 
energy of expansion. At the same time the strong shock 
which marks the boundary of the remnant converts this 
macroscopic kinetic energy of expansion back into micro- 
scopic internal form. Thus there is a continuous recycling 
of the original explosion energy between the micro and 
macro scales. This continues until either the external pres- 
sure is no longer negligible compared to the internal, or 
the time-scales become so long that radiative cooling be- 
comes important. The time scales for the conversion of ki- 
netic energy to internal energy and vice versa are roughly 
equal and of order the dynamical time scale of the rem- 
nant which is of order the age of the remnant, hence the 
approximately self-similar evolution. 

In terms of particle acceleration the theory assumes 
that strong collisionlcss shocks in a tenuous plasma auto- 
matically and inevitably generate an approximately power 
law distribution of accelerated particles which connects 
smoothly to the shock-heated particle distribution at 'the- 
rmal' energies and extends up to a maximum energy con- 
strained by the shock size, speed, age and magnetic field. 
The acceleration mechanism is a variant of Fermi accelera- 
tion based on scattering from magnetic field structures on 
both sides of the shock. A key point is that these scatter- 
ing structures are not those responsible for general scat- 
tering on the ISM, but strongly amplified local structures 
generated in a non-linear bootstrap process by the acceler- 
ated particles themselves. As long as the shock is strong it 
will be associated with strong magnetic turbulence which 
drives the effective local diffusion coefficient down to val- 
ues close to the Bohm value. As pointed out by Achterberg 



et al. (1994) the extreme sharpness of the radio rims of 
some shell type SNRs can be interpreted as observational 
evidence for this type of effect. The source of free energy 
for the wave excitation is of course the strong gradient 
in the energetic particle distribution at the edge. Thus in 
the interior of the remnant, where the gradients are absent 
or much weaker, we do not expect such low values of the 
diffusion coefficient. 

The net effect is that the edge of the remnant, as far 
as the accelerated particles are concerned, is both a self- 
generated diffusion barrier and a source of freshly accel- 
erated particles. Except at the very highest energies the 
particles produced at the shock are convected with the 
post-shock flow into the interior of the remnant and effec- 
tively trapped there until the shock weakens to the point 
where the self-generated wave field around the shock can 
no longer be sustained. At this point the diffusion barrier 
collapses and the trapped particle population is free to 
diffuse out into the general ISM. 

In terms of bulk energetics, the total energy of the ac- 
celerated particle population is low during the first ballis- 
tic phase of the expansion (because little of the explosion 
energy has been processed through the shocks) but rises 
rapidly as t « t$w During the sedov-like phase the to- 
tal energy in accelerated particles is roughly constant at 
a significant fraction of the explosion energy (0.1 to 0.5 
typically). However, this is because of the dynamic recy- 
cling described above. Any individual particle is subject 
to adiabatic losses on the dynamical time-scale of the rem- 
nant, while the energy lost this way goes into driving the 
shock and thus generating new particles, distributed over 
the whole energy spectrum. 

3. Spallation reactions within SNRs 

3.1. Qualitative overview 

We now turn to the production of Li, Be and B (LiBcB) by 
spallation reactions within a SNR. As emphasized above, 
there are two obvious mechanisms. One is the irradiation 
of the CNO ejecta by accelerated protons and alphas. It 
is clear that the fresh CNO nuclei produced by the SN 
will, for the lifetime of the SNR, be exposed to a flux of 
energetic particles (EPs) very much higher than the aver- 
age interstellar flux, and this must lead to some spallation 
production of light elements. This process starts at about 
^sw with a very intense radiation field and continues with 
an intensity decreasing roughly as i?~ 3 oc t~ 6 / 5 (where R 
is the radius of the SNR) until the remnant dies. 

The second process is that some of the CNO nuclei 
from the ejecta are accelerated, either by the reverse shock 
in its brief powerful phase at t ~ t$w or by some of 
this material managing to get ahead of the forward shock. 
This later possibility is not impossible, but seems unlikely 
to be as important as acceleration by the reverse shock. 
Calculations of the Raleigh- Taylor instability of the con- 
tact discontinuity do suggest that some fast-moving blobs 
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of ejecta can punch through the forward shock at about 
£sw(Jun & Norman 1996), and in addition Ramaty and 
coworkers have suggested that fast moving dust grains 
could condense in the ejecta at t < tsw an d then pen- 
etrate through into the region ahead of the main shock. 
In all these pictures acceleration of CNO nuclei takes place 
only at about tsw and the energy deposited in these accel- 
erated particles is certainly less than the explosion energy 
-EsNj although it might optimistically reach some signif- 
icant fraction of that value (say 10%). Crucially the 
accelerated CNO nuclei are then confined to the interior 
of the SNR and will thus be adiabatically cooled on a 
rather rapid timc-scalc, initially of order tsw- 

3.2. Evaluation of the first process (forward shock) 

From the above arguments, it is clear that SNe do induce 
some Be production. Now the question is : how much? Let 
us first consider the irradiation of the ejecta by particles 
(H and He nuclei) accelerated at the forward shock during 
the Sedov-like phase - process 1. We have already indi- 
cated that detailed studies of acceleration in SNRs show 
that the fraction of the explosion energy given to the EPs 
is roughly constant during the Sedov-like phase and of or- 
der 0.1 to 0.5 or so. Let d\ be that fraction. Since the EPs 
are distributed more or less uniformly throughout the in- 
terior of the remnant, the energy density can be estimated 
as 



£cr 



(1) 



where R is the radius of the remnant and -Esn is the ex- 
plosion energy. 

To derive a spallation rate from this we need to assume 
some form for the spectrum of the accelerated particles. 
Shock acceleration suggests that the distribution function 
should be close to the test-particle form f(p) oc p~ 4 and 
extend from an injection momentum close to 'thermal' 
values to a cut-off momentum at about p max = 10 5 GeV/c. 
The spallation rate per target CNO atom to produce a Be 
atom is then obtained by integrating the cross sections 



'-'spall 



■)f(p)dp, 



(2) 



Pth 



with the normalisation J E(p)f(p)4:Trp 2 dp = £cr- Look- 
ing at graphs of the spallation cross-sections for Be (as 
given, e.g., in Ramaty et al. 1997), it is clear that these 
cross-sections can be well approximated as zero below a 
threshold at about 30-40 MeV/n and a constant value 
Co — 5 x 10 -27 cm 2 above it. One then obtains roughly : 



'spall 



(JO 

^CR -, 

mc 1 



1 - In (p th /mc) 
ln(pmax/mc) 



0.2a o c^, (3) 
mc 



where pth — mc/5 is the momentum corresponding to the 
spallation threshold and m refers to the proton mass. For- 
tunately, for this form of the spectrum the upper cut-off 



and the spallation threshold only enter logarithmically. A 
softer spectrum would lead to higher spallation yields and 
a stronger dependence on the spallation threshold. 

Using Eq. (|l|) and the adiabatic expansion law for the 
forward shock radius, R = i?sw(*Asw) 2 ^ 5 , we can now 
estimate the total fraction of the CNO nuclei which will 
be converted to Be during the Sedov-like phase as : 
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QiEsn Po , 



mc 



( tsw 

\ ^ond 



1/5 



(4) 
(5) 

(6) 



where as usual po denotes the density of the ambient 
medium into which the SNR is expanding and M e j is 
the total mass of the SNR ejecta. We now recall that the 
sweep-up time is given in terms of the SN parameters and 
the ambient number density, no ~ Po/ni, as 



-1/3 / o h \ 1/3 

n ( _3_ Mq x 

4-7T m 



(7) 



where t> C j ~ (2EsN/M e j) 1 ^ 2 is the velocity of the ejecta, 
or numerically : 



t S w = (1.410 3 yr) 



/ M oj V/ E s 



\10MqJ V10 51 erg 



V 1cm -3 / 



(8) 



Replacing in Eq. (||) and using canonical values of 
E SN = 10 51 erg and Af cj = 10 Af Q , we finally get : 



no \ 2 / 3 



^4xl0- lo M-^4 
V 1 cm a 



\ t en d ) 



1/5 



(9) 



Clearly this falls short of the value of order 10 -8 re- 
quired to explain the observations, even for values of 0\ 
as high as 0.5. It might seem from Eq. (^) that very high 
ambient densities could help to make the spallation yields 
closer to the needed value. This is however not the case. 
First, the above estimate does not take energy losses into 
account, while both ionisation and adiabatic losses act to 
lower the genuine production rates. Second, and more sig- 
nificantly, the ratio ^swAend (and a fortiori its fifth root) 
becomes very close to 1 in dense environments, lowering 0i 
quite notably (see Fig. ^|). In fact, it turns out that there 
is no Sedov-like phase at all in media with densities of or- 
der 10 4 cm -3 , the physical reason being that the radiative 
losses then act on a much shorter time-scale, eventually 
shorter than the sweep-up time. 
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3.3. Evaluation of the second process (reverse shock) where U = |7ri? 3 e is the total kinetic energy of the par- 

. • i , ■ • • • / 1 2 _ • . i ■ . • • . • 



Let us now turn to the second process, namely the spal- 
lation of energetic CNO nuclei accelerated at the reverse 
shock from the SN ejecta and interacting within the SNR 
with swept-up ambient material. We have argued above 
that this reverse shock acceleration is only plausible at 
times around t$w and certainly the amount of energy 
transferred to CNO nuclei cannot be more than a frac- 
tion of Esn ■ Let #2 be the fraction of the explosion energy 
that goes into accelerating the ejecta at or around tsw, 
and #cno the fraction of that energy that is indeed trans- 
ferred to CNO nuclei. These particles are then confined to 
the interior of the remnant where they undergo spallation 
reactions as well as adiabatic losses. Let us again assume 
that the spectrum is of the form f(p) oc p~ 4 . Then the 
production rate of Be atoms per unit volume is approxi- 
mately 



0.2 



npo-Q £cr 
mc 14 



(10) 



where £qr now refers to the accelerated CNO nuclei, 
the factor 14 comes from the mean number of nucleons 
per CNO nucleus and the factor 0.2, as before, from the 
f(p) oc p~ 4 spectral shape (assuming the same upper cut- 
off position, but this only enters logarithmically). Integrat- 
ing over the remnant volume, we obtain the spallation rate 
at tsw '■ 



cLV B . 
dt 



0.2cr c 



#CNO#2-EsN 

14mc 2 



(11) 



Now the adiabatic losses need to be evaluated rather 
carefully. It is generally argued that they act so that the 
momentum of the particles scales as the inverse of the 
linear dimensions of the volume occupied. Accordingly, in 
the expanding spherical SNR the EPs should lose momen- 
tum at a rate p/p = —R/R, reminiscent, incidentally, of 
the way photons behave in the expanding universe. In our 
case, however, the situation is complicated by the fact that 
the EPs do not push directly against the 'walls' limiting 
the volume of confinement, which move at the expansion 
velocity, V = R, but are reflected off the diffusion barrier 
consisting of magnetic waves and turbulence at rest with 
respect to the downstream flow, and thus expanding at 
velocity 

To see how this influences the actual adiabatic loss 
rate, it is safer to go back to basic physical laws. Adia- 
batic losses must arise because the EPs are more or less 
isotropised within the SNR and therefore participate to 
the pressure. Now this pressure, P, works positively while 
the remnant expands, implying an energy loss rate equal 
to the power contributed, given by : 



AU 
~dt 



F v = 



PdS x-R= -3ttR 2 RP, 



(12) 



tides. Considering that P = |e in the non-relativistic 
limit (NR) and P = |e in the ultra-relativistic limit (UR), 
Eq. ( |l2"| ) can be re-writen as : 



de 
dt 



-M < UR >- 



(13) 



Finally, dividing both sides by the space density of the 
EPs and noting that E — p 2 /2m in the NR limit, and 
E = pc in the UR limit, we obtain the momentum loss 
rate for individual particles, valid in any velocity range : 



p 
P 



3R 
AR' 



(14) 



From this one deduces that at the time when the rem- 
nant has expanded to radius R, only those particles whose 
initial momenta at tsw were more than (i?/i?sw) 3 ^ 4 Pth 
are still above the spallation threshold. For a p~ A distri- 
bution function the integral number spectrum decreases 
as p~ x and thus the number of accelerated nuclei still ca- 
pable of spallation reactions decreases as i? -3 / 4 oc i~ 3 / 10 . 
For a softer accelerated spectrum the effect would be even 
stronger because there are proportionally fewer particles 
at high initial momenta. 

This being established, we can integrate Eq. (|ll|) over 
time, to obtain the total production of Be atoms : 



A/Be = 0.2n<7o c 



#CNO#2-EsN 

14mc 2 



end / t 

\*sw 



-3/10 



that is 



?CNO#2-EsN 

14mc 2 



/ £cnd\ 

\tswj 



7/10 



dt (15) 



(16) 



Dividing by the total number of CNO nuclei in the 
ejecta, A oj;C NO = (6 , CNo/14)A cji tot ^ (6» CNO A4j/14m), 
we get the final result : 



#2-Esn Po , 
oqc r7T f SW 



uic 



2 / tend\ 

7 V*sw/ 



7/10 



1 - 



^end / 



7/10 



(17) 



Note that we assumed that the mass fraction of CNO in 
the ejecta is the same as the energy fraction of CNO in 
the EPs (which was the original meaning of #cno)- Con- 
sidering that all nuclear species have the same spectrum in 
MeV/n, and thus a total energy proportional to their mass 
number, this simply means that the acceleration process is 
not chemically selective, in the sense that the composition 
of the EPs is just the same as that of the material passing 
through the shock. 
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Numerically, again with E'sn = 10 erg and M t 



10 Mq, we finally obtain 

b 2 ~ i x io- w e 2 1 n ° 



1 cm 



2/3 



Vsw/ 



7/10 



1 - 



tsw 



7/10' 



(18) 



3-4- Relative contribution of the two processes 

It is worth emphasizing the similarity between expressions 
(^) and (|l7]) that we obtained for the spallation rates per 
CNO nuclei by the two processes considered here. This 
formal analogy allows us to write down their relative con- 
tributions straightforwardly : 



h &2 2 7 / 10 

— = -Q- X ~ (tend tsw) — 



1 - (^SwAcnd) 7 ^ 10 



1 



(isw/^end) 1 ^ 5 



(19) 



As is often the case, this similarity is not fortuitous 
and has a physical meaning. The two processes may indeed 
be regarded as 'dual' processes, the first consisting of the 
irradiation of the SN ejecta by the ambient medium, and 
the second of the ambient medium by the SN ejecta. The 
'symmetry' is only broken by the dynamical aspect of the 
processes. First, of course, the energy imparted to the EPs 
in both cases needs not be the same, for it depends on 
the acceleration efficiency as well as the total energy of 
the shock involved (forward or reverse). This is expressed 
by the expected ratio 62/ Ox- And secondly, in the first 
process one has to fight against the dilution of the ejecta 
- integration of (t/tsw)~ 6 ^ 5 , see Eq. (||) - while in the 
second process one fights against the adiabatic losses - 
integration of (t/tsw)~ 3 ^ W , see Eq. (|l5|). This is expressed 
by the last factor in Eq. (|l^). 

Clearly the latter decrease of the production rates is 
the least dramatic, and the reverse shock process must 
dominate the LiBeB production in supernova remnants. 
However, this conclusion still depends on the genuine effi- 
ciency of reverse shock acceleration, and once the relative 
acceleration efficiency 82/81 is given, the weight of the first 
process relative to the second still depends on the total 
duration of the Sedov-like phase, appearing numerically 
in Eq. ( [i"9| ) through the ratio WdAsw, which in turn de- 
pends on the ambient density, no- The expression of tsw as 
a function of the parameters has been given in Eq. (@), so 
we are left with the evaluation of the time, t nd, when the 
magnetic turbulence collapses and the EPs leave the SNR, 
putting an end to Be production. We argued above that 
i cn d should correspond to the end the Sedov-like phase, 
when the shock induced by the SN explosion becomes ra- 
diative, that is when the cooling time of the post-shock 
gas becomes of the same order as the dynamical time. 

In principle, the cooling rate can be derived from the 
so-called cooling function, A(T)(erg cm 3 s _1 ), which de- 




Fig. 1. Cooling function (bold line) as a function of the 
temperature for a medium with metallicity lower than 
~ 10~ 2 Z Q . The dashed lines illustrate the graphical de- 
termination of t cn d, for an ambient density uq = 10 cm -3 
(see text). 

pends on the physical properties of the post-shock mate- 
rial, notably on its temperature, T, and metallicity, Z : 



nA(T) ' 



(20) 



where n is the post-shock density, equal to 4no if the com- 
pression ratio is that of an ideal strong shock (nonlinear 
effects probably act to increase the compression ratio to 
values larger than 4). As for the dynamical time, we sim- 
ply write 



^ R 

Tdyn ~ 



5 

-t. 

2 



(21) 



To obtain t e nd> we then need to solve the following 
equation in the variable t, obtained by equating r coo i and 
r dyn given above : 

3k B T 



t 



20n A(T) ' 



(22) 



where it should be clear that the right hand side also 
depends on time through the temperature, T, and thus 
indirectly through the cooling function too. In the non- 
radiative SNR expansion phase, the function T(t) is ob- 
tained directly from the hydrodynamical jump conditions 
at the shock discontinuity : 



T 



3m 



-V 1 



(23) 



or numerically 



T « (2 x 10 5 K) 



Ei 



SN 



10 51 



erg 



lem 



if t 



\ 10 5 yr 



(24) 



To solve Eq. (p2|), we still need to know the cooling 
function A(T). In the range of temperatures correspond- 
ing to the end of the Sedov-like phase, 10 5 K < T < 10 7 K, 
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it happens to depend significantly on metallicity, with dif- 
ferences up to two orders of magnitude for metallicities 
from Z = to Z = 2Z (Bohringer & Hensler 1989). Be- 
cause we focus on Be production in the early Galaxy we 
adopt the cooling function corresponding to zero metallic- 
ity, represented in Fig. |l| (adapted from Bohringer & Hens- 
ler 1989), which holds for values of Z up to ~ 10~ 2 Z Q . 

For high enough ambient densities, the shock will be- 
come radiative early in the SNR evolution, when the tem- 
perature is still very high, say above T Jb 2 10 6 K. In this 
case, the cooling function is dominated by Bremsstrahlung 
emission and can be written analytically as : 

A Br (T) a (2.4 10- 23 erg cmV 1 ) f-J^ J . (25) 

Substituting from (§J) and @ in Eq. (§|) and solving 
for t, we find : 

To check the consistency of our assumption A « Ae r 
(i.e. T Z 210 6 K), let us now report Eq. (§|) in © 
and write down the temperature T^ n d at the end of the 
Sedov-like phase : 

T End a(2xlO^K)(-|^) 1/4 ( T ^) 1/2 , (27) 
\10 oi erg/ Vlcm / 

which means that the above analytical treatment is valid 
only for ambient densities greater than about 100 cm -3 . 
For lower densities, we must solve Eq. p2] ) graphically. 
First, we invert Eq. ( p4| ) to express t as a function of tem- 
perature, then we plot the function /(T) = ik-gT / (2Qn^t) 
on the same graph as A (see Fig. [j] for an example), find 
the value of T at intersection, and finally convert this value 
into the sought time i cn d making use again of Eq. (|24|). The 
results, showing i cn d as a function of the ambient density, 
are shown in Fig. ^. 

We now have all the ingredients to plot the efficiency 
ratio of the two processes calculated above. Figure || shows 
the ratio </>2/</>i given in Eq. ( [i"9| ) as a function of the am- 
bient density, assuming that 9\ = 9 2 . Two different values 
of the ejected mass have been used, corresponding to dif- 
ferent progenitor masses (~ 10 — 40 M Q ). It can be seen 
that low densities are more favourable to the reverse shock 
acceleration process. This is due to iondAsw being larger, 
implying a larger dilution of the ejecta (process 1 less effi- 
cient) and smaller adiabatic losses, which indeed decrease 
as t^ 1 (process 2 more efficient). The part of the plot 
corresponding to 4>2/4>i < 1 is not physical, because it re- 
quires i on d < isw, which simply means that the Sedov-like 
phase no longer exists and the whole calculation becomes 
groundless. Note however that in Fig. [3] the energy im- 
parted to the EPs has been assumed equal for both pro- 
cesses, which is most certainly not the case. Actually, if 
2 /9 1 = 0.1 (e.g. 9 1 = 10% and 6>i = 1%), then process 1 



10 6 f 




10~ 1 10° 1 1 1 2 1 3 1 4 1 s 10 s 

ambient density (cm" 3 ) 

Fig. 2. Comparison of tend and isw as a function of the 
ambient density, no, for two different values of the mass 
ejected by the supernova (10 and 30 M Q ). The dashed line 
shows the asymptotic analytic estimate of Eq. ( p6| ) . 
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Fig. 3. Comparison of the Be (and B) production effi- 
ciency through the forward and reverse shock accelera- 
tion processes, for two values of the ejected mass (10 and 
30 Mq). The ratio (see text) is plotted as a func- 

tion of the ambient density, assuming that both processes 
impart the same total energy to the EPs (0\ =62)- 

is found to dominate Be production during the Sedov-like 
phase, regardless of the ambient density. 

4. Spallation reactions after the Sedov-like phase 

At the end of the Sedov-like phase, the EPs are no longer 
confined and leave the SNR to diffuse across the Galaxy. 
At the stage of chemical evolution we are considering here, 
there are no or few metals in the interstellar medium 
(ISM), so that energetic protons and a particles acceler- 
ated at the forward shock will not produce any significant 
amount of Be after i cn d (although Li production will still 
be going on through a+a reactions). In the case of the sec- 
ond process, however, the EPs contain CNO nuclei which 
just cannot avoid being spalled while interacting with the 
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(H+He)/(C+0) <0.1 




10 -2 10-' 10° 10 1 

1 H / 4 He source abundance ratio 



Fig. 4. Production efficiency of Be, as a function of the 
EP composition (all abundances are by number). The or- 
dinate is the number of Be nuclei produced by spallation 
reactions per erg injected in the form of EPs. A thick tar- 
get has been assumed, with zero metallicity. Carbon and 
Oxygen abundances were set equal in the EP composition. 

ambient H and He nuclei at rest in the Galaxy. This may 
be regarded as a third process for Be production, which 
lasts until either the EPs are slowed down by Coulombian 
interactions to subnuclear energies (i.e. below the spalla- 
tion thresholds) or they simply diffuse out of the Galaxy. 
Since the confinement time of cosmic rays in the early 
Galaxy is virtually unknown, we shall assume here that 
the Galaxy acts as a thick target for the EPs leaving the 
SNR, an assumption which actually provides us with an 
upper limit on the spallation yields. 

Unlike the first two processes evaluated above, this 
third process is essentially independent of dynamics. Thus, 
time-dependent calculations are no longer needed and, 
from this stage on, the calculations made by Ramaty et al. 
(1997) or any steady-state calculation is perfectly valid. In 
particular, the ambient density has no influence on light 
element production, since a greater number of reactions 
per second, as would result from a greater density, im- 
plies an equal increase of both the spallation rates and the 
energy loss rate. Once integrated over time, both effects 
cancel out exactly, and in fact, given the energy spectrum 
of the EPs, the efficiency of Be production (and Li, and 
B), expressed as the number of nuclei produced per erg in- 
jected in the form of EPs, depends only on their chemical 
composition. 

Results are shown in Fig. || for different values of the 
source abondance ratios, H/He and (H + He)/(C + O), al- 
lowing one to derive the spallation efficiency for any com- 
position. Two-steps processes (such as 12 C -I- 1 H — » 10 B 
followed by 10 B + 1 H — > 9 Be) have been taken into ac- 
count. Test runs show good agreement with the results of 
Ramaty et al. (1997). 

As can be seen on Fig. |I| pure Carbon and Oxygen 
have a production efficiency of about 0.22 nuclei/erg, while 
this efficiency decreases by at least a factor of 10 for com- 




3 V— r — — — — — n 

-icr 1 io° io 1 10 2 10 3 10 4 io 5 
ambient density (cm A -3) 

Fig. 5. Fraction of the energy imparted to the EPs at time 
tsw which is still available at i e nd) after the Sedov-like 
phase, as a function of othe ambient density. 

positions with hundred times more H and He than metals 
(or about ten times more by mass) . According to models of 
explosions for SN with low metallicity progenitors, the av- 
erage (H+He)/(C + 0) ratio among the EPs should indeed 
be expected to be ^ 200, unless selective acceleration oc- 
curs to enhance the abundance of the metals. As a conse- 
quence, efficiencies greater than ~ 10 -2 nuclei/erg should 
not be expected, so that a production of ~ 4 1 48 atoms 
of Be requires an energy of ~ 4 1 50 erg to be imparted to 
the EPs. This seems very unlikely considering that the to- 
tal energy available in the reverse shock (the source of the 
EPs) should be of order one tenth of the SN explosion en- 
ergy, not to mention the acceleration efficiency. Moreover, 
a significant fraction of the energy originally imparted to 
the EPs has been lost during the Sedov-like phase of the 
SNR evolution through adiabatic losses. 

To evaluate the 'surviving' fraction of energy, it suffi- 
cies to go back to Eq. ([l4]), which indicates that when the 
radius of the shock is multiplied by a factor rj, the momen- 
tum p of all the particles is multiplied by a factor r]~ 3 / 4 . 
It is worthwhile noting that, because of their specific mo- 
mentum dependence, adiabatic losses do not modify the 
shape of the EP energy spectrum. In our case, f(p) ex p -4 , 
so that when all momenta p are divided by a factor £, the 
distribution function f(p) is divided by the same factor £. 
To see that, the easiest way is to work out the number of 
particules between momenta p and p + dp after the mo- 
mentum scaling. This number writes dN' = f (p)4wp 2 dp, 
where f'(p) is the new distribution function. Now diV' 
must be equal to the number of particles that had mo- 
mentum between Qp and £(p + dp), which is, by definition, 
dN = /(Cp)47r(Cp) 2 Cdp. Equating dN and dN' yields the 
result f'(p)=(-\f(p). 

Putting all pieces together, we find that when the 
shock radius R is multiplied by a factor r), the distribu- 
tion function and, thus, the total energy of the EPs are 
multiplied by ?7~ 3 / 4 . Now considering that R increases as 
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i 2 / 5 during the Sedov-like phase, we find that the total 
energy of the EPs decreases as £~ 3 / 10 . Note that this is 
nothing but an other way to work out the decrease of the 
spallation rates for our second process during the Sedov- 



like phase (cf. Sect. 3.3). Finally, we find that a fraction 
(iendAsw)" 3 ^ 10 of the initial energy imparted to the EPs 
is still available for spallation at the end of the Sedov-like 
phase. This factor is plotted on Fig. || as a function of 
the ambient density. It can be seen that for no = 1 cm -3 , 
the energy available to power our third process of light 
element production has been reduced by adiabatic losses 
to not more than one third of its initial value, and less 
than one half for densities up to 100 cm" 3 . Clearly, high 
densities are favoured (energetically) because they tend to 
shorten the Sedov-like phase, and therefore merely avoid 
the adiabatic losses. 

5. Discussion 

Since light element production in the interstellar medium 
obviously requires a lot of energy in the form of supernu- 
clear particles (i.e. with energies above the nuclear thresh- 
olds) as well as metals (especially C and O), it is quite 
natural to consider SNe as possible sources of the LiBcB 
observed in halo stars. We have analysed in detail the spal- 
lation nucleosynthesis induced by a SN explosion on the 
basis of known physics and theoretical results relating to 
particle shock acceleration. Two major processes can be 
identified, depending on whether the ISM or the ejecta 
are accelerated, respectively at the forward and reverse 
shocks. In the first case, the EPs consist mostly of pro- 
tons and alpha particles and must therefore interact with 
C and O nuclei, which are much more numerous within 
the SNR than in the surrounding medium (especially at 
early stages of Galactic evolution). The process will thus 
last as long as the EPs stay confined in the SNR, i.e. ap- 
proximately during the Sedov-like phase, but not more. 
In the second case, freshly synthesized CNO nuclei are ac- 
celerated, and Be production occurs through interaction 
with ambient H and He nuclei. The process is then di- 
vided into two, one stretching over the Sedov-like phase, 
with the particles suffering adiabatic losses, and the other 
one occuring outside the remnant, with only Coulombian 
losses playing a role. 

We have calculated the total Be production in these 
three processes, taking the dynamics of the SNR evolu- 
tion into account (dilution of the ejecta by metal-poor 
material and adiabatic losses). The results are shown in 
Fig. H for processes 1 and 2 (from Eqs. (||) and (|l8|)). We 
find that with canonical values of 9\ = 0.1, 62 = 0.01, 
M C j = 10 Mq and a mean ambient density no = 10cm -3 , 
the fraction of freshly synthesized CNO nuclei spalled into 
Be in these processes is </>i ~ 810~ n and </> 2 ~ 310~ n , 
respectively, which is very much less than the value 're- 
quired' by the observations, discussed in the introduction 
(0obs ~ 3 10 -8 ). Even allowing for unreasonably high val- 
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Fig. 6. Spallation efficiency of CNO during the Sedov-like 
phase, as a function of the ambient density. The fraction of 
freshly synthesized CNO nuclei being spalled to Be by pro- 
cesses 1 (4>i) and 2 (^> 2 ) is obtained from Eqs. (||) and (18) 
and the values of i on d derived in Sect. 3.4, for two values 
of the ejected mass (10 and 30 Mq). 



ues of the acceleration efficiency, Q\ ~ 62 ^> 1, the total Be 
production by processes 1 and 2 would still be more than 
one order of magnitude below the observed value. 

As suggested by Fig. ||and our analytical study, higher 
densities improve the situation. However, even with no — 
10 3 cm~ 3 and acceleration efficiencies equal to 1, the Be 
yield is still unsufficient. Moreover, it should be noted that 
our calculations did not consider Coulombian energy losses 
(because they are negligeable as compared to adiabatic 
losses for usual densities), which become important as the 
density increases and therefore make the Be yield smaller. 
Finally, since we are trying to account for the mean abun- 
dance of Be in halo stars, as compared to Fe, we have to 
evaluate the Be production for an ambient density corre- 
sponding to the mean density encountered around explo- 
sion sites in the early Galaxy, which is very unlikely to 
be as high as 10 3 cm~ 3 . It could even be argued that al- 
though the gas density might have been higher in the past 
than it is now (hence our 'canonical value' no = 10 cm" 3 ), 
the actual mean density about SN explosion sites could be 
lower than 1 cm -3 , because most SNe may explode within 
superbubble interiors, where the density is much less than 
in the mean ISM. 

Thus, our conclusion is that processes 1 and 2 both 
fail in accounting for the Be observed in metal-poor stars 
in the halo of our Galaxy. Concerning the third process, 
adopting canonical values for the parameters again leads 
to unsufficient Be production, as noted in the previous 
section. While higher densities improve the situation by 
avoiding the adiabatic energy losses, one should neverthe- 
less expect at least half of the EP energy to be lost in this 
way, for any reasonable density. This means that even if 
10% of the explosion energy is imparted to EPs acceler- 
ated at the reverse shock, which is certainly a generous 
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upper limit, the required number of ~ 4 1 48 nuclei of Be 
per SN implies a spallation efficiency of ~ 0.1 nucleus/erg. 
Now Fig. |^ shows that this requires an EP composition in 
which at least one particle out of ten is a CNO nucleus. In 
other words, the ejected mass of CNO must be of the order 
of that of H and He together. None of the SN explosion 
models published so far can reproduce such a requirement, 
and so there is clearly a problem with Be production in 
the early Galaxy. 

The results presented here are in fact interesting in 
many regards. First, they show that it is definitely very 
difficult to account for the amount of Be found in halo 
stars. Consequently we feel that the main problem to 
be addressed in this field of research is probably not the 
chemical evolution of Be (and Li and B) in the Galaxy, as 
given by the ensemble of the data points in the abundance 
vs metallicity diagrams (e.g. whether Be is proportional to 
Fe or to its square) but, to begin with, the position of any 
of these points. Are we able to describe in some detail one 
process which could explain the amount of Be (relative to 
Fe) present in any of the stars in which it is observed? 
The answer, we are afraid, seems to be no at this stage. 
It is however instructive to ask why the processes investi- 
gated here have failed. Concerning process 1 (acceleration 
of ISM, interaction with fresh CNO within the SNR), the 
main reason is that the CNO rich ejecta are 'too much 
diluted' by the swept-up material as the SNR expands, 
so that the spallation efficiency is too low (or the avail- 
able energy is too small). However, it seems rather hard 
to think of any region in the Galaxy where the concen- 
tration in CNO is higher than inside a SNR during the 
Sedov-like phase (especially in the first stages of chemical 
evolution)! So the conclusion that process 1 cannot work, 
even with a 100% acceleration efficiency, seems to rule out 
any other process based on the acceleration of the ISM, 
initially devoided of metals. 

The other solution is then of course to accelerate CNO 
nuclei themselves, which provides the maximum possible 
spallation efficiency, independently of the ambient metal- 
licity. Every energetic CNO will lead to the production 
of as much Be as possible given the spallation cross sec- 
tions and the energy loss rates. The latter cannot physi- 
cally be smaller than the Coulombian loss rate in a neu- 
tral medium, and this leads to the efficiency plotted in 
Fig. U Unfortunately, a significant amount of the CNO 
rich ejecta of an isolated SN can only be accelerated at the 
reverse shock at a time around the sweep- up time, tsw- 
This means that i) the total amount of energy available 
is smaller than the explosion energy (probably of order 
10%, i.e. ~ 10 50 erg), and ii) the accelerated nuclei will 
suffer adiabatic losses during the Sedov-like phase, reduc- 
ing their energy by a factor of 2 or 3. As shown above, 
this makes process 2-3 incapable of producing enough Be, 
as long as the EPs have a composition reflecting that of 
the SN ejecta. 



This suggest that a solution to the problem could be 
that the reverse shock accelerates preferentially CNO nu- 
clei rather than H and He. For example, recent calcula- 
tions have shown that such a selective acceleration arises 
naturally if the metals are mostly condensed in grains (El- 
lison et al. 1997). The proposition by Ramaty et al. (1997) 
that grains condense in the ejecta before being accelerated 
could then help to increase the abundance of CNO in the 
EPs. However, we have to keep in mind that any selective 
process called upon must be very efficient indeed, since as 
we indicated above, the data require that the EP compo- 
sition be as rich as one CNO nuclei out of ten EPs, which 
means that CNO nuclei must be accelerated at least ten 
times more efficiently than H and He. This would have 
to be increased by another factor of ten if the energy ini- 
tially imparted to the EPs by the acceleration process were 
only a factor 2 or 3 lower (i.e. ~ 3 1 49 erg, which is more 
reasonable from the point of view of particle acceleration 
theory). Clearly, more work is needed in this field before 
one can safely invoke a solution in terms of selective ac- 
celeration. 

As can be seen, playing with the composition to in- 
crease the spallation efficiency has its own limits, and in 
any case, Fig. ^ gives an unescapable upper limit, obtained 
with pure Carbon and Oxygen (at least for the canonical 
spectrum considered here - other spectra were also investi- 
gated, as in Ramaty et al. (1997), leaving the main conclu- 
sions unchanged). This would then suggest that another 
source of energy should be sought. However, the constraint 
that it should be more energetic than SNe is rather strong. 

Another interesting line of investigations could be the 
study of the collective effects of SNe. Most of the mas- 
sive stars and SN progenitors are believed to be born 
(and indeed observed, Melnik & Efremov 1995) in asso- 
ciations, and their joint explosions lead to the formation 
of superbubbles which may provide a very favourable en- 
vironment for particle acceleration (Bykov & Fleishman 
1992). Parizot et al. (1998) have proposed that these su- 
perbubbles could be the source of most of the CNO-rich 
EPs, and Parizot & Knoedlseder (1998) further investi- 
gated the gamma-ray lines induced by such an energetic 
component. The most interesting features of a scenario 
in which Bc-producing EPs are accelerated in superbub- 
bles is that i) when a new SN explodes, the CNO nuclei 
ejected by the previous SNe are accelerated at the for- 
ward shock, instead of the reverse shock in the case of an 
isolated SN, which implies a greater energy, and ii) no 
significant adiabatic losses occur, because of the dimen- 
sions and low expansion velocity of the superbubble. This 
makes the superbubble scenario very appealing, and it will 
be investigated in detail in a forthcoming paper. 

However that may be, we should also keep in mind 
that when we say that a process does not produce enough 
Be, it always means that it does not produce enough Be 
as compared to Fe. Now it could also be that SN explo- 
sion models actually produce too much Fe. The point is 
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that Be is compared to Fe in the observations, while it 
has no direct physical link with it. Indeed, Be is not made 
out of Fc, but of C and O. So to be really conclusive, the 
studies of spallative nucleosynthesis should compare the- 
oretical Be/O yields to the corresponding abundance ra- 
tio in metal-poor stars. Unfortunately, the data are much 
more patchy for Be as a function of [O/H] than as a func- 
tion of [Fe/H], especially in very low metallicity stars. The 
usually assumed proportionality between O and Fe could 
turn out to be only approximate, as recent observational 
works possibly indicate (Israelian et al. 1998; Boesgaard 
et al. 1998; these observations, however, still ?? need to be 
confirmed by an independent method, all the more that 
they come into conflict with several theoretical and obser- 
vational results; cf. Vangioni-Flam et al. 1998b). We shall 
address this question in greater detail in the attending 
paper (Parizot and Drury, 1999, Paper II). 

Finally, we wish to stress that the calculations pre- 
sented in this paper rely on a careful account of the dy- 
namics of the problem. More generally, time-dependent 
calculations are required to properly evaluate the spalla- 
tion processes in environments where compositions and 
energy densities are evolving. In particular, as argued in 
Parizot (1998), no variation with density can be obtained 
with a stationary model, since an increase in the density 
induces an equivalent and cancelling increase in the spal- 
lation rates and the energy loss rates. By contrast, we 
have shown that all three of the processes considered here 
are more efficient at higher density - a result which could 
not have been found otherwise. Detailed, numerical time- 
dependent calculations will be presented in paper II, with 
conclusions similar to those demonstrated here. 
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